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Abstract. In this paper we consider a scalar transport equation with constant 
coefficients on domains with discrete space and continuous, discrete or general time. 
We show that on all these underlying domains solutions of the transport equation 
can conserve sign and integrals both in time and space. Detailed analysis reveals 
that, under some initial conditions, the solutions correspond to counting stochastic 
processes and related probability distributions. Consequently, the transport equation 
could generate various modifications of these processes and distributions and provide 
some insights into corresponding convergence questions. Possible applications are 
suggested and discussed. 
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1. Introduction 

Scalar transport equation with constant coefficients u t + ku x = belongs among the simplest 
partial differential equations. Its importance is based on the following facts. Firstly, it describes 
advective transport of fluids, as well as one-way wave propagation. Secondly, it serves as a base for 
a study of hyperbolic partial differential equations (and is consequently analysed also in numerical 
analysis). Thirdly, its nonlinear modifications model complex transport of fluids, heat or mass. 
Finally, its study is closely connected to conservation laws (see [5] or [T3"]). 

Properties and solutions of partial difference equations have been studied mainly from numerical 
(e.g. [13]) but also from analytical point of view (e.g. [5]). Meanwhile, in one dimension, there 
has been a wide interest in the problems with mixed timing, which has been recently clustered 
around the time scales calculus and the so-called dynamic equations (see |3], [H]). Nevertheless, 
there is only limited literature on partial equations on time scales (see [T] , [3] , [IB] ) ■ These papers 
indicate the complexity of such settings and the necessity to analyze basic problems like transport 
equation. Our analysis is also closely related to numerical semidiscrete methods (e.g. [T3] Section 
10.4]) or analytical Rothe method (e.g. [17]). 

In this paper we consider a transport equation on domains with discrete space and general 
(continuous, discrete and time scale) time (see Figure [T]). We show that the solutions of transport 
equation does not propagate along characteristics lines as in the classical case and feature behavior 
close to the classical diffusion equation. Our analysis of sign and integral conservation discloses 
interesting relationship between the solutions on such domains and probability distributions related 
to Poisson and Bernoulli stochastic processes. These counting processes are used to model waiting 
times for occurence of certain events (defects, phone calls, customers' arrivals, etc.), see [2], [TU] or 
15J for more details. Consequently, considering domains with general time, we are able not only 
to generalize these standard processes but also generate transitional processes of Poisson-Bcrnoulli 
type and corresponding distributions. Moreover, our analysis provides a different perspective on 
some numerical questions (numerical diffusion) and relate it to analytical problems (relationship 
between the CFL condition and regressivity) . Finally, it also establishes relationship between 
the time scales calculus and heterogeneous and mixed probability distributions in the probability 
theory. 

In Section [3] we summarize well-known features of the classical transport equation. In Section 
[4]we consider a transport equation with discrete space and continuous time. In Section[5]we solve 
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Figure 1. Examples of various domains considered in this paper. We study 
domains with discrete space and continuous (Section [4]), discrete (Section [5| and 
general time (Section [6j). 



the problem on domains with discrete time. In Section [6] we generalize those results to domains 
with a general time and prove the necessary and sufficient conditions which ensure that the sign 
and both time and space integrals are conserved (Theorem [l8| ) . Finally, in Section [7] we discuss 
convergence issues, applications to probability distributions and stochastic processes and provide 
two examples. 



2. Preliminaries and Notation 

The sets R, Z, N denote real, integer and natural numbers. Furthermore, let us introduce 
N = N U {0} and = [0, oo). Finally, we use multiples of discrete number sets, e.g. a-multiplc 
of integers is denoted by aZ and defined by aZ = {. . . , —2a, —a, 0, a, 2a, . . .}. 

Partial derivatives are denoted by Ut(x,t) and u x {x,t) and partial differences by 



(1) 



A t u(x,t) 



u(x, t + /if) — u(x, t) 



Hi 



and W x u{x,i) 



u(x,t) - u(x - (M x ,t) 



where \i t and \x x denote step sizes in time and space. 

In Section [6j we consider time to be a general time scale T, i.e. an arbitrary closed subset 
of K. Time step could be variable, described by a graininess function /j,t ■ T —> Rq . We use 
the partial delta derivative u At which reduces to Ut in points in which fitit) = or to A t u in 
those t in which fj,(t) > 0. Similarly, we work with the so-called delta-integral which corresponds 
to standard integration if T = R or to summation if T = Z. Finally, the dynamic exponential 
function e p (x,Xo) is defined as a solution of the initial value problem (under the regressivity 
condition 1 +p(i)/i(i) ^ 0) 

ix A (t)=p(t)x(t), 
\x(0) = l. 

For more details concerning the time scale calculus we refer the interested reader to the survey 
monograph [3]. 

Given function u(x,t), by u(x, •) we mean functions of one variable having the form u(0,t), 
it(l, t), etc. Similarly, by w(-, t) we understand one-dimensional sections of u(x, t) having the form 
u(x, 0), u(x, 1), etc. 
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(2) 



3. Classical Transport Equation 

Let us briefly summarize essential properties of the classical transport equation 

Ju t (x,t) + ku x (x,t) = 0, (6R + ,ieR, 
1 u(x,0) = (j)(x), xeR. 

with (f> E C 1 . Typical features whose counterparts are studied in this paper include: 

• the unique solution u(x, t) — <j)(x — kt) could be obtained via the method of characteristics, 
the solution is constant on the characteristic lines where x — kt = C, 

• consequently, the solution conserves sign, i.e. if 4>{x) > then u(x, t) > 0, 

oo 

• moreover, the solution conserves integral in space sections, i.e. if J <f>{x)x — K, then 



/ 



u(x, f)x = K, for all t > 0, 



• finally, the solution conserves integral in time sections in the following sense. For k > 
we have that 

OO X 

u(x, t)t — — / <p(s)s. 
k J 

— oo 

Consequently, if 4>(x) = for x > xq, then the integral along time sections is constant for 
all x > Xq. 

4. Discrete Space and Continuous Time 

In contrast to the classical problem ^ we consider a domain with discrete space and the 
problem 

{ut{x,t) + kV x u{x,t) = 0, t G Kq", x £Z, 
U(X ' 0) = l0, x^O, 
where A > 0, k > and V x u reduces tcfj] 

Vjiifa;, i) = u(x, t) — u(x — 1, t). 
One could rewrite the equation in ([3| into 

Ut(x, t) = —ku(x, t) + ku(x — l,t), 
which implies that the problem ([3| could be viewed as an infinite system of differential equations, 
lem 1. The unique solution of the problem ^ has the form: 



(4) u(x,t) 



A * t x e -M^ t G Mj,:z; G N , 



x 



0, t G R$,x g Z,x < 0. 



Proof. First, let us observe that u(x,t) — for all t G K.q , x < 0. The uniqueness of the trivial 
solution for x < follows e.g. from [151 Corollary 1] or more generally from [51 Theorem 3.1.3]. Let 
us prove the rest (i.e. x > by mathematical induction. Obviously, we have that u(0,t) — Ae~ kt , 
since u t (0, t) = -ku(0, t) + ku(-l, t) = -ku(0, t) and u(0, 0) = A. 



^We assume that k > so that the solution is bounded and does not vanish. Moreover, we use the nabla 
difference instead of delta difference. The single reason is the simpler form of the solution J3J. If we used the delta 
difference, we would consider k < and the solution would propagate to the quadrant with t > and x < 0. This 
applies also to the problems which we study in the following sections. 
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• 


u(x,l) 


■ 


u(x,4) 


▲ 


uU.8) 


♦ 


ulx.lll 



FIGURE 2 . Solution of the transport equation with discrete space and continuous 
time ([3]) with A — 1 and k = 1. 



□ 



Moreover, if we assume that neNo and u(x,t) = A^jt x e kt , then u(x + l,t) satisfies 

J u t (x + 1, t) = -ku(x + 1, t) + A^Pe"' 1 *, fel + , 
|u(a;,0) = 0. 

One could use the variation of parameters to show that the unique solution is u(x + 1, i) 
A-j^Yy,£ :E+1 e~ fct , which proves the inductive step and consequently finishes the proof. 

Let us analyze the sign and integral preservation of ([3|. 

lem 2. The solution of the problem ^ conserves the sign, the integral in time and the sum in 
the space variable. 

Proof. The sign preservation follows from the positivity of all terms in Q. Next, we could use 
integration by parts to obtain (we skip the details since we prove this result in more general 



settings in Theorem 15 ) 
(5) 

Similarly, summing over x we get 



u(x, t)t 



(6) 



e4< 



x e -kt 



Ae - Kt = Ae- kt e kt = A. 



A 
k' 

(kt) x 



□ 

If we go deeper and analyze values obtained in ^ and ^ we get the first indication of the 
relationship of the semidiscrete transport equation with stochastic processes. 

Remark 3. If A = k then time sections of the solution Q generate the probability density function 
of Erlang distributions ( note that for x — we get the exponential distribution and that Erlang 
distributions are special cases of Gamma distributions). 

Similarly, if A = 1 the space sections of Q form the probability mass functions of Poisson 
distributions. 

Consequently, if A = k = 1 the solution u(x, t) describes Poisson process. All these facts are 
further discussed in Section [?| 

We conclude this section with two natural extensions. Firstly, we mention possible generaliza- 
tions to other discrete space structures. 
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Remark 4. If we consider the problem ^ on a domain with a discrete space having the constant 
step fi x > 0, not necessarily fi x — 1, we obtain qualitatively equivalent problem, since 

Ut[x, t) + k = Ut(X, t) H \u(x, t) — u[x — fi x ,t)) 

fix \l>x 

— Ut{x, t) + k\7 x u(x, t). 

In contrast to the rest of this paper, the value of fi x does not play essential role here. Therefore, 
for presentation purposes, we restricted our attention to fi x = 1 . 

Finally, we discuss more general initial condition and show that the solution is the sum of point 
initial conditions which justifies their use not only in this section but also in the remainder of this 
paper. 



Corollary 5. The unique solution of 

ju t (x,t) + kV x u(x,t) = 0, t eR^,x eZ, 
\u(x,0) = C x . 



(7) 

is given by 

(8) »(*,*)= £ C t P^e^. 

^-^ [x — i)\ 

i— — oo ' 

Proof. One could split ^ into problems with point initial conditions, use Lemma [TJ to solve them 
and then employ linearity of the equation to get □ 

5. Discrete Space and Discrete Time 

In this section, we assume that both time and space are homogenously discrete with steps 
fit > and fi x > respectively. In other words, we consider a discrete domain 

Q = {(x,t) = (mfi x , nfi t ) , with m G Z, n G No} . 

The transport equation and the corresponding problem then have the form 

{A t u(x,t) + W x u(x,t) = 0, (x,t)eSl, 
( ^ l A > x = > 
U(X ' 0) = l0, x^O, 

where A > 0, k > 0. Using the definition of partial differences in |l]), we can easily rewrite the 
equation in ^ into 

(10) u(x, t + fi t ) = (l - — ] u{x, t) + — u(x - fi x ,t) 

\ fix J fix 

and derive the unique solution. 

lem 6. Let m G Z and n G No- The unique solution of ^ has the form: 



At" " ^\ , n>m>0, 



(11) u(mfi x ,nfi t ) = 



l^x J \ l^x 

0, < n < m, or m < 0. 

Proof. First, let us show that the solution vanishes uniquely for u(—mfi x , nfi x ) = 0, for all m, n G 



N. Consulting ( 10 1, we observe that the value of u{—mfi x , nfi x ) is obtained as a linear combination 
of initial conditions u(—mfi Xl 0), u{— (m+ l)fi x , 0), u{— (m+n)/x x , 0), i.e. a linear combination 
of n + 1 zeros. 

We prove the rest of the statement by induction. Apparently, 

u (0,n/xt)=fl-— V«(0,0)=a(i k ^ 
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Next, let us assume that u(mfi x ,nfi t ) satisfies (111, then 

knt\ n 



u((m + l)fi x ,nfi t ) = I 1 ) u((m + l)/x x ,0) 

r m+ i=0 v TO/ v TO/ v "' r m =0 r 2 =0ri=0 

/ , \ n—(m+l) /, \ m+1 n—l r 3 — 1 r 2 — 1 

= a(i-^) (^) E •••EE 1- 

" ' 7 r m+1 =0 r 2 =0r 1= 



At this stage, let us observe the properties of the falling factorials (see e.g. [H Section 2.1] or 
[T2l Section 2.1]) to get that 

n-1 rs-lra-X f n \ 

^ " ' ^ ^ 1 ~ (m + 1)! ~~ \m + 1 J ' 

r m+ i=0 r 2 =0ri=0 V / \ / 

which finishes the proof. □ 

The closed-form solution enables us to analyze sign and integral conservation, 
lem 7. If the inequality 

(Dl) 1 - > 0, 

ftoZds then the solution of ^ satisfies 

(i) u(x,t) > 0, 

oo 

(ii) E u(m/j, x ,t) is constant for all t = {0, /it, ■ ■}, 

m— — oo 

oo 

(iii) E w(x, n/zt) is constant for all x = {0, fi x , 2fi x , . . .}. 

71=0 

Proof. (i) The inequality follows immediately from Lemma [6] . 
(ii) If we fix t and sum up the equation (|10[) over 2 we get 



^ u(mfi x , t + fit) = ( 1 -) u(mfi x ,t)-\ u((m- l)fi x ,t). 

\ (^X J f^X 

ra= — oo x 7 m— — oo m=— oo 

The assumption (Dl) implies that the sum on the left hand side is a linear combination of 
two sums on the right hand side. Since these sums are equal, we get that 

oo oo 

^ u(mfi x ,t + fit) = ^ u(mfi x ,t). 

m=—oc m= — oo 



(iii) Similarly, one could sum up the equation ( |10[ ) over t to get for a fixed x > 

OO / ^ \ OO oo 

rz./z t ) =11 ) u(s, n/i t ) H u(a; - /Xa, nfi t ). 



n—l v / n—0 n—0 

Since u(x, 0) = for x > we have that 

oo oo 

^ u(x, nfit) = ^u(x - fb X) nfi t ). 

?i=0 n=0 

Once again, we could study the solutions' relationship to probability distributions. 



□ 



thm 8. Let u(x,t) be a solution of Then the space and time sections fi x u(x, •) and fi t u(-,t) 
form probability mass functions if and only if the assumptions (Dl), 
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(D3) 
hold. 



A\i x = 1, 



Proof. Lemma [7] yields that the solutions are nonnegative and conserve sums. It suffices to include 
step-lengths fi x and /i t and identify conditions under which fi x J2 X u(x, 0) = 1 and [it J2t *) = 
1. Given the initial condition, the former sum is equal to Afi x . Hence the assumption (D3). Finally, 

since u(—fi x ,t) = 0, the equation (10 1 implies that u(0,nfi t ) = A (l — l ^ ± ^j ■ Consequently, 



1 = Aii 



OO 

E 



i 



\i x 



A\i x 



□ 



Corollary 9. Let u(x,t) be a solution of (j9J). Then the space and time sections fi x u(x,-) and 
Htu(-, t) form probability mass functions if and only if k = 1, fi t < [i x and A = 

Proof. (D2) and (D3) hold if and only if k — 1. Consequently, (Dl) could be satisfied if and only 
if fi t < Ms- D 

Closer examination again reveals that the sections form probability mass functions of discrete 
probability distributions related to Bernoulli counting processes. 



Remark 10. Let us consider the solution (111. If we put A = k = fi x = 1 and \x t = p we get 

u(n, m • p) = ( J (1 — p) n m p m , n > m. 

which forms, for each fixed n 6 No, a probability mass function of the binomial distribution. 
Similarly, for each fixed m £ No, p — fit-multiple forms a probability mass function of a version of 
the negative binomial distribution (the value p-u(n,m-p) describes a probability that for m failures 
we need n trials). Consequently, the solution of ([9| describes a counting Bernoulli stochastic 
process (see ^\). 
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6. Discrete Space and General Time 

Let us extend the results from the last two sections by considering more general time structures. 
Let T be a time scale such that minT = and supT = +oo. In this paragraph we consider domains 



and the problem: 



Q = {(x,t) : x e fi x Z,te T} . 



u At (x,t) + kV x u(x,t) =0, (x,t)eQ, 



where A > 0, k > and V x u(x, t) is the backward difference defined in and u At is the delta- 



derivative in time variable. Since the space is discrete, we could again rewrite the equation in 12 
into 

k 

(13) u At (x, t) = (u(x, t) — u(x — fJ> x ,t)) . 

In order to conserve the sign of solutions we assume that 



(TS1) 1 - > 0i 

f^x 

i.e. the condition which is similar to the positive regressivity in the time scale theory (e.g. [H 
Section 2.2]) or the so-called CFL condition in the discretization of the transport equation (e.g. 
P31 Section 4.4] 



Let u be a solution of (12). One could use Q31 Proposition 5.2] to show that u(x,t) — for all 



x < is the unique solution there. Since u(—fx x ,t) — 0, we could see that u At (0,t) = — — it(0, i). 



Given the initial condition and assumption (TS1), we get u(0,t) = Ae_ _k_(t;Q), where e_j^(t;0) 



is a time scale exponential function (see [H Section 2]). 



lem 11. The solution of (12) satisfies 
(i) lim u(0, f) = 0, 



t— too 



(ii) Ju(0,t)At = A^. 



Proof. (i) Follows directly from the assumption (TS1) and the properties of the exponential 
functions [H Section 2.2], 

(ii) 



/ u(0,t)At = A e_ ft (t;0)At = lim -A^- (e_j^(t;0) - l) = A 

J J f»x *->00 k V A»x / 



k 



□ 



Unique solutions of u{mpL x ,t) could be found using the variation of constants (see e.g. [H 
Theorem 2.77]). However, these computations depend critically on a particular time scale and 
cannot be performed in general. For example, one could compute that the second branch of the 
solution has the form 



k At 



u(fj, x ,t) = A — e_j^(t;0) 



Hx Jo 1- 



tlx 

This implies that we can't derive closed-form solutions as in previous sections. Formally, these 
solutions can be expressed as Taylor-like series with generalized polynomials whose form depends 
on particular time scale (see [14] and [H Section 1.6]). We determine these solutions in special cases 
(see Lemmata [lj [uj and |l9|) . Therefore, we are forced use another means to show the properties of 
solutions we are interested in. 
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lem 12. Let x G /i^N. // (TS1 ) is satisfied and u(x — fi x , t) > and for allt € T and u(x — fj, xi t) > 
at least for one t € T then u{x, t) > /or ai! t e T. 



Proof. First, note that u(x,0) = for all x > 0. Consequently, (13 1 implies that u At (x,t) > at 
the beginning of the support of u{x — fi x ,t) and u(x,t) is strictly increasing there. 

• If t is right-scattered then we can rewrite the equation ( 13 ) into 



u(x,t + fJL t ) = I 1 - 



U(x, t) + — u(x - fi x ,t). 



If u(x,t) > 0, then this is the weighted average of two nonnegative values and thus 
nonncgative as well. 
• If t is right-dense then the equation ( 13 1 has the form 

u t (x, t) = u{x, t) H u(x — fi x ,t). 

t^x \^x 

Since both u(x — ^ x ,t) > and u(x,t) > 0, we have that Ut(x,t) > —-^-u(x,t) and thus 
u(x, t) cannot become negative. 

Following the induction principle (e.g. [H Theorem 1.7]), we could see that u(x,t) > for all 



t G T. 



□ 



Lemma [12] serves as the inductive step in the proof of the sign-conservation. 



thm 13. // ( |TST| ) holds then u{x,t) > for all (x,t) e Cl. 

Proof. We prove the statement by mathematical induction. Firstly, u(0,t) = Ae__^(t;0) > 0. 

Mas 

Secondly, if u(x, t) > then Lemma 12 implies that u(x + fi x ,t) > which finishes the proof. □ 

The following auxiliary lemma shows that the variation of constant formula which generates 
further branches of solutions conserve zero-limits at infinity. 

lem 14. Let us consider a time scale T, a constant K such that 1 — fiK > and a function 

oo 

/ : T — > [0, oo) such that the integral J f(t)At is finite. If we define g : T — > [0, oo) by 

o 



g(t)= / e. K (t,a(r))f(r)AT, 



then lim g{t) = 0. 



Proof. Since J f(t)A(t) is finite we know that for each e > there exists T > such that for all 
o 

t e T, t > T the inequality 



(14) 



OO 

J /(r)Ar < 



holds. Similarly, properties of time scale exponential function imply that for each e > and T > 
there exists R> T such that for all £ € T, t > R the following inequality is satisfied 



(15) 



e_ K (t;a{T))AT < 



2F' 
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with F = max /(f). Consequently, inequalities (14) and (15) imply that for each for each e > 

there exists T > and R > T such that for alii > R 

t 



g(t)= / e_ K (t,a(T))f(T)A T 



T t 

e^ K (t;a(T))f(T)AT+ [ e_ K (t; <r(r))/(r)Ar 



1 t 

<Fj e_ K (t;a(r))Ar + J /(r)Ar 



< F 



2F 2 



which implies that lim g(t) = 0. 

t— >oo 

Consequently, we are able to show that the integrals are constant for each fixed x > 0. 



□ 



thm 15. If (TS1) holds and u(x,t) is a solution of (12) then 

oo 

u(x,t)At = J u(0,t)At = A- 



for all x G yU^No- 

Proof. We proceed by mathematical induction. 



• For x — the convergence of the integral to 4y follows from Lemma 11 (ii). 

• Let us fix x £ fi x N and assume that the statement holds for a function u(x — /i x ,t). If we 



(16) 



integrate ([13 ) we get 

u At (x, t)At 



k 



u(x, t)At 



u(x - (i x ,t)At 



Let us concentrate on the left-hand side term. The variation of constants formula ([H 
Theorem 2.77]) implies that 

/ 

u(x,t) — / e__k_ (i, <j{t))u(x — /i x , t)At. 



Consequently, Lemma 



14 



implies that lim u(x,t) — 0. Using the initial condition u(x, 0) 

t— f oo 



0, we could rewrite theTcft-hand side of ( 16 ) into 



u '(i,t)At= lim u(x, t) — u(x, 0) = 0. 

t— >oo 



This implies that ( 16 1 could be rewritten into 
k 



= 



or equivalently into 



^ u(x,t)At-J u(x - ^,t)At| , 



u(x, t)At 



u(x - (i x ,T)Ar, 



which finishes the proof. 
Finally, we show that the integrals (sums in this case) remains constant in time as well. 



□ 
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/ X x >5 
U •< 



• 


u(x,l) 


■ 


u(x,4) 


▲ 


u(x.8) 


♦ 


ulx.lll 



FIGURE 4. Solution of the transport equation with discrete space and general 



time Q with A = 1, /j, x = 1, k = 1 and T = [j [i, i + |] . 

i=0 



thm 16. If (TS1) holds and u{x,t) is a solution of (12 1 then 



J oo 

/ u(x, t)Ax = Hx) j u(mfi x , t) = A[i x 



/or oi/ i G T. 



Proof. Let tt(a;, t) be a solution of (12). We define a function 5 : T — > R by 



^ oo 

:= / i)Ax = ^ u(mfi x ,t), 

J ™ — n 



and show that S A *(t) = for all t e T. 
We can rewrite the equation in (fl2| into 



Consequently, 
(17) 



u At (x, t) = u(x, t) H u(x — (j, x , t). 

H"x l^x 



s At (t) = i± x ^ uAt ' *) 



m=0 
oo 



(18) 
(19) 



= -k ^ v>(m(j, x ,t) + k ^ u (( m ~ i)/^, t) 



m=0 



m=0 



= 0. 



We have to justify the first equality (17 1, i.e. the interchangability of the delta-derivative and 



summation at each to G T. If to is right-scattered the non-negativity of the solution implies 

oo oo 

u(miJ, x ,to + fi t (to)) - J2 u{mfi x ,t a ) 



S At (t ) 



m=0 



m=0 



Mt(*o) 

Eu(mfi x ,t + fi t (to)) -u(mn x ,to) f a„ . \ 

m=0 m=0 



Mt(*o) 

If to is right-dense and there is a continuous interval [to,s], s > to, we show that the sum 



"^2 u * (m\x x , t) converge uniformly on [to, s]. First, let us note that ( 18 ) yields that this is implied 

m=0 
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by the uniform convergence of ^2 u(mfj, x ,t). One could use Corollary 5 to get (k = k//j, x ): 

m=0 

m=0 m=0 \ i=0 ^ / 

= e -«(t-t ) " t o)) m . c . 

i— 1 ml ^ *' 



If X)i=o ^ s finite (i.e. 5(to) is finite), then this sum converge uniformly on an arbitrary closed 
interval. Finally, if t is right-dense and there is no continuous interval [to, s], s > to, we consider a 
function v(x, t) with v(m/j, x , to) = u(mfj, x , to) for all m such that v is a solution on a domain with 
a continuous interval [to,s], s > to- Obviously, the equation in (12) implies that Vt(mfj, x ,to) = 
u At (miJ, x ,to) for all m. Moreover for each 8 > there is 9 > such that for all t € [to, to + @}t'- 

oo oo oo 

(1 - S) >^ v(mfi x ,t) < 2J u{mp x ,t) < (1 +£) v(m{i Xl t). 

rn—0 m—0 m— 

Consequently, 

CXD CXD 

= X! u At (mfi x ,t ) = 22 v t {mfi x ,t ) 

m—0 m—0 
/ oo \ / oo \ A * 



Taking into account the fact that u(x, 0) is given by the initial condition in (12), we see that 
5(0) = Afi x . Consequently, (|l7|)-(fl9| imply that S(t) = Afj, x . □ 

We could now study the relationship with probability distributions and we begin by generalizing 
probability density and mass functions. We say that a function / : T — > Kq is a dynamic 
probability density function if 



/oo 
f(t)At = 1. 
-oo 



Note that if T = K then / is a probability density function. If T = /i t Z then fi t f is a probability 
mass functions (see Therorem [8]) . 

Combining Theorems [15] and [16] we get the necessary and sufficient condition for sections to 
generate probability distributions. 



lem 17. Let u{x,t) be a solution of (12) 

r>rn co-fit i ftiin ^rtnn nil t* d 1 1 T^J„ t~f nn A r\v\ In if - 

k 



(1) u(-,t) is a dynamic probability density function for all x £ Mz^o if an d only if = 1 



and fJ>t(t) < [i, x for all t 6 T. 
(2) u(x, ■) is a dynamic probability density function for all t 6 T if and only if Afj, x = 1 and 
( TSlj ) holds. 

Proof. The proof is a direct application of Theorems [15] and [l6[ □ 
Finally, we provide the necessary and sufficient condition for both sections. 



thm 18. Let u{x,t) be a solution of (12). Then both u(x, •) and u(-,t) are dynamic probability 

density functions for all t € T and x £ if and only if k — 1, Afi x = 1 and m(t) < fi x for 
each t € T. 

Proof. The proof follows from Lemma |17| □ 
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U(;t) 


u(0,-) 


u(x, ■), x > 


Zxl 


Poisson dist. 


exponential dist. 


Erlang (Gamma) dist. 


ZxpZ 


binomial dist. 


geometric dist. 


negative binomial dist. 



Table 1 . Correspondence of time and space sections with probability distributions 



7. Applications 

As suggested in Remarks [3] and [lO] the time and space sections of solutions of the transport 
equation on various domains generate important probability distributions (cf. Table [lj. 

In other words, the solutions correspond to the so-called counting stochastic processes describing- 
number of occurrences of certain random events (arrival of customers in a queue, device failures, 
phone calls, scored goals, etc.) (e.g. [TH Chapters 4 and 5], [10]). They have following properties 

(1) probability of number of events (occurrences) at time t is given by u(-,t) (Poisson distri- 
bution, binomial distribution), 

(2) probability distribution of the time of the first occurrence is given by u(0, t) (exponential 
or geometric distribution), 

(3) probability distributions that at least x events have happened until time t are given by 
u(x — 1, •) (Erlang or negative binomial distributions), 

(4) probability distribution of the waiting time until the next occurrence is given by u(0, t) 
(exponential or geometric distribution). 

Our analysis in Section [6] summarized in Theorem [Ls] suggests that properties (l)-(3) are 
conserved on general domains Z x T. Properties (2)-(3) are conserved in the sense of Definition 



?? (see Examples 20 and [21] below) . Property (4) does not apply because of the underlying 
inhomogeneous time structure. 

The convergence relationship between the distributions from Table [I] is well-known [10] . Our 
analysis strengthens this relationship since the convergence is based on the solution of the same 
partial equation with changing underlying structures. 

We conclude this section by suggesting two applications which emphasize the time scale choice. 
First, let us consider Bernoulli trials with non-constant probability of successes. For example, [7] 
shows that the probability that a goal is scored in each minute of the association football match is 
not constant but increases throughout the game, especially in the last minutes of each half-time. 
Let us derive an explicit solution on arbitrary heterogeneous discrete structure. 

lem 19. Let us consider a heterogenous discrete time scale T = {0, + [L2-, • ■ • , S"=i IHi ■ ■ •}■ 



Then the solution of (12 1 has the form 



(20) u £ m = A ]T II K i iL i~" 




where Ki = 1 — k^-, Li = kj^- and Pf! denote a set of all permutation vectors containing q ones 
and r zeros. 

Proof. We base our proof on the relationship 

u I m^ x ,y^/Xj ] = ( 1 - u{x,t) + f^u(x - (i x ,t) 

and proceed by induction. First, the initial condition implies that the statement holds for n = 0. 



Next, let us assume that the statement holds for n € No, i.e. (|20|) is satisfied. Then we have 

u I 

n+l 



yn(j, x , YH=i f^nj = for m ^ {0, 1, . . . , n + 1}. Furthermore, for m = 
f n+l \ 

0, IM = K n+1 (AKiK 2 ...K n ) + = AK,K 2 . . . K n K n+1 . 



u 

\ i=l 
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Next, for me (1, 2, . . . , n): 

(n+l \ n n 

mn x ,J2m)= Kn+iA II K i iL l~ ni + L n + iA £ ]J K?L\~*< 

t=i / neP%r m 4=1 7I "eP , ;i™ +1 i=1 

n+l 

Finally, for m = n + 1 we have 

/ n+l \ 

u I (n + 1)^, 2J Mi = + £ n+ i (ALiL 2 . . . L n ) = AL Y L 2 . . ■ L n L n+1 . 



□ 



We could immediately apply this result to obtain generalizations of standard Bernoulli pro- 
cesses. 

exmp 20. Heterogeneous Bernoulli process. Let us consider a repeated sequence of trials 
and assume that the probability of success pi in i-th trial is non-constant, in contrast to standard 
Bernoulli process discussed in Section^ If we construct a discrete time scale 



T = < 0,pi,pi +p 2 , ■ ■ -,^2pi 



i=l 



n \ n /i\ 7r ;/-\ 1_7r i 



i/ien t/ie solution u(x,t) of (12 1 generates the probability distributions discussed above. Let us 
choose, for example, A = p, x = k = 1. Then, u(-, X^—i Pi) *s f/ie probability mass function 
describing number of successes in the first n trials. Moreover, piu{x, •) is the probability mass 
function of the number of trials needed to get x + 1 successes. 

To illustrate, let us choose T — /o, §j 5 + §>••• j X] i+T' ' • "J ^° s ^ u( ki a process in which the 

probability of successful trial decreases harmonically. We could use Lemma \T^ to determine that: 

1— Wi 

x — ^ \ x — t — r / ^ \ / * ^ 

U I m, 

For example, the probability mass function for the first successful trial appearing in k-th trial, i.e. 
Piu(0,-), has the form 

™-wh)> fceN - 

Finally, we consider a mixed time scale, which, coupled with the transport equation, generates 
mixed processes and distributions. 

exmp 21. Stop-Start Bernoulli-Poisson Process. Let us assume that a device is regularly 
used throughout a constant period and then switched off for another one. Let us assume that the 
probability of failure when the device is in use is determined by a continuous process, whereas the 
probability of failure in the rest mode is given by a discrete process. This leads to mixed probability 
distributions which could be generated e.g. by 



T=U 



i=0 



1, 1 + 



Again u{x, ■) describes the mixed probability distribution ofx + 1 failures, in the sense of Definition 
11. Similarly, u(-,t) is the probability mass function describing the number of failures at time t. 
Note that the probability of failure in the rest mode is given by the length of the discrete gap (cf. 
Definition 11). As in the previous example, we are not able to find the closed-form solutions but 
one could tediously solve the separate equations to get that: 

«(<M)= ^e 5 ^' 
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, N 2t + n a , 
u(l,t) = - e* , 

. , At 2 + Ant + (n 2 - An) „ , 



2! • 2 n + 2 

3! • 2™+ 3 



. , 8t 3 + I2nt 2 + 6(n 2 - An)t + (n 3 - I2n 2 + 16n) 



, . polynomial of order x » , 

u(a;,t) = — r-r e"*, 



/or ?i G No (n-th continuous part) and t £ [ra, n + 5] . See Figure^] /or illustration. 

8. Conclusion and Future Directions 

There is a number of open questions related to the anaysis presented in this paper. In Section 
[6] we were unable to provide a general closed-form solution of problem (12). With the connection 



to probability distributions, is it possible to provide one for further special choices of T (see e.g. 
Examples [20] and [21])? 

In the classical case, the solution is propagated along characteristics. Obviously, our analysis in 
Sections |4] and [5] implies that this is not the case on semidiscrete domains. However, one could show 
that at least the maxima are propagated along characteristics on discrete-continuous or discrete- 
discrete domains (computing directly or using modes of probability distributions). Having no 
closed- form solutions on time scales, could we prove this property for an arbitrary time scale? 
This question is closely related to modes of the corresponding probability distributions and the 
question could be therefore formulated in more general way. Can we, at least in special cases, 
determine the descriptive statistics related to the generated probability distributions? 

From the theoretical point of view, there is also a natural extension to consider a transport 
equation with continuous space and general time, or general space and time. The applicability 
of this settings is limited by the fact that such problems does not conserve sign in general (cf. 



assumption u t (t) < u x in Theorem 18) 
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